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Abstract 

We study the Zee-Babu twodoop neutrino mass generation model and look for a possible 
flavor symmetry behind the tri-bimaximal neutrino mixing. We hnd that there probably exists 
the ji-r symmetry in the case of the normal neutrino mass hierarchy, whereas there may not be 
in the inverted hierarchy case. We also propose a specific model based on a Froggatt-Nielsendike 
Z§ symmetry to naturally accomplish the \i-t symmetry on the neutrino mass matrix for the 
normal hierarchy case. 
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I. INTRODUCTION 



Neutrino oscillation experiments have almost completely established that neutrinos 
have tiny masses and mix with each other through the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) leptonic mixing matrix l|. From the latest global analysis of three-neutrino 
mixing 2|, one currently has the following best fit values with la errors: 



Am 2 21 



(7.59 ±0.20) x 10" 5 eV 2 , 

—2.36 ± 0.11) x 10 -3 eV 2 for inverted hierarchy 
+2. 46 ±0.12) x 10" 3 eV 2 for normal hierarchy 
sin 2 12 = 0.319 ± 0.016, sin 2 9 23 = 0.462+° 



Am li 



)+0.082 
1.050' 



sin 2 0i3 



0.0095i°;°J? 



The data indicate the existence of, at least, two massive neutrinos with a very suggestive 
neutrino mixing matrix, that is, the tri-bimaximal (TB) mixing matrix j^]: 

1 2 y/2 N 
v tb = ^= -1 V2 ■ (2) 

\ 1 -V2 V3j 

However, the standard model (SM) neither includes neutrino mass terms nor provides us 
with any explanation for the TB mixing. Clearly, we need new physics beyond the SM. 
' n fact, many extensions of the SM have been proposed so far. For instance, in the type-I 
4], type-II [5] and type-Ill j^] seesaw mechanisms, the SM is extended by introducing 
extra heavy fermions or scalars to generate neutrino masses suppressed by the mass scale 
of the heavy particles, while in Ref. (3] tiny neutrino masses come from the dimension- 
five Weinberg operators. These scenarios have been extensively studied with some flavor 
symmetries to explain the TB mixing 

Yet another possibility of leading to tiny neutrino masses is to use radiative corrections. 
It was first pointed out by Zee in Ref. 13[ in which new scalars are added in the Higgs 
sector with neutrino masses induced at the one- loop level. After that, a two- loop scenario 
called Zee-Babu model [lij] was proposed 1 . In these kinds of scenarios, discussions about 
the neutrino phenomena can be much different form those of the tree level scenarios 



1 Other types of multi-loop scenarios have also been studied in Ref. [1 
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because the induced neutrino mass matrix elements are the products of the anti-symmetric 
Yukawa coupling and charged lepton mass. Hence, it is non-trivial whether a flavor 
symmetry can play an important role in the radiative scenarios. Since there is a claim 
that the original Zee model may not be able to reproduce current neutrino oscillation data 



16[, we focus on the Zee-Babu two-loop model in this Letter. We re-analyze the model 



and try to explain the TB pattern of neutrino mixings in terms of the \i-t symmetry 
which is the prime candidate of a flavor symmetry in the tree level scenarios. Note that 
other phenomenological studies have been discussed in Refs. 171. Il8l|. 

This Letter is organized as follows. In Section II, we summarize the Zee-Babu model 
and show some definitions of parameters. In Section III, we investigate the model along 
with the yU-r symmetry. We propose a specific flavor model in Section IV. Finally, we 
conclude our discussions in Section V. 



II. ZEE-BABU MODEL 

In addition to the SM particles, the Zee-Babu model contains two SU(2)l singlet new 
scalars: a singly charged scalar h^ 1 and doubly charged scalar /c ±=t . Accordingly, new 
interactions appear and terms relevant to our study are 



C ZB = F ab (L a CL b h + ) + Y ab (£l a C£ m k ++ ) - fih + h + k" + h.c 



(3) 



where C is the charge conjugation matrix, L a=e ^^ T stand for the left-handed SU(2)l 
doublet leptons and are the right-handed singlet charged leptons in the diagonal basis 
of the charged lepton mass matrix. F ab and Y ab are 3x3 complex Yukawa matrices, 
parametrized as 



ab 







fefi fi 

U o i 

\ -fer -fur J 



1/J.T 



Y„i 



^ Dee Ve^i Ver 

Ue^l y^L^L UjlT 
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(4) 



The Majorana neutrino mass term is induced at the two-loop level as depicted in Fig. [T] 
with the mass matrix given by 



Mab = 8n(fa C m c y* d m d f db )I, 



(5) 
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FIG. 1: Two-loop diagram for Majorana neutrino masses. 



where m a indicate the charged lepton masses and 

1 1 W'^/ T,'!:/ 1 ,, , , .og ( 6) 



■ (16ir 2 ) 2 Ml Jo Jo *x+[(M t /Mtf-l]y + y 2 6 x + (Mt/Mtfy 
is the two-loop integral function with the masses of the new scalar s, Mk and M^. Note 
that Eq. ([6]) is simplified by neglecting the charged lepton masses [171 ] . The elements of 
the neutrino mass matrix in Eq. (JHJ) are written as 



M 33 = %nfl T {-u^ - 2f eT u eil - ft T u ee )I, 



(7) 



-^12 = 8/J,f* T (-f eT U TT - /e^W^T + fe^fer^er + fe^e^I = M 2 1, 
M n = 8/J,f^ T (f eT U^ T + ] e ^^ + /eVer + fejer^e^l = M 3 1, 
M 2 3 = 8flf^ T (uj^ T + fer^er ~ fe^ey. ~ fejer^ee)! = M 32 , 

with the following redefinitions of parameters: 

fen = 7^ Jer = -J— , w ab = m a y ab m b . (8) 

It is clear that the mass matrix in Eq. (jSJ) always has a zero-eigenvalue because of 
the vanishing determinant of F^. Although all three active neutrinos should have non- 
zero masses if we take into account the higher order loop contributions, we ignore these 
contributions in this Letter. 

Furthermore, we can embed the terms associated with the electron mass into u TT , upl- 



and 00^ terms, such that 



UJ 



/'/' 



Then, we obtain the following simplified mass matrix 

M ab = 8fiflu'^IM ab 

with 



(9) 



(10) 



M, 



ah 



Cfl 



fer^rr feuU^ur fer^nr ft 



c l- 1 



\ 



* 



UJ 



)1T 



where Un T and uj tt are defined as 



UJ 



ijj 



fJT 



OJ-r 



ijj 7 



UJ 



UJ 



(12) 



U 



PMNS 



ci 3 si 3 e" 
1 

!<5 



^ I Ci 2 S12 0\ / 1 (A 



/ 



-Sl2 C12 

Ij 



e^' 2 

V° i; 



(13) 



As partially discussed in Ref. [17|, we can represent f efl , feri^nr an d cu TT in terms of 
the neutrino mass ratios, mixing angles and CP violating phases in the PMNS matrix, 
parametrized by 

d o o W 

C23 S23 
V -s 23 c 23 J \-s 13 e 16 ci 3 
where 5 and 7 are the Dirac and Majorana CP phase, respectively, and Sij(cij) = 
sin 6ij (cos dij) > 0. Since we consider the diagonal basis of the charged leptons, the Ma- 
jorana mass matrix is diagonalized by the PMNS matrix, such that Up MNS M. Upmns — 
diag(mi, m% tus). In the case of the normal mass hierarchy, the four parameters are 
described as 



fefj, 


Sl2 S 2 3 

Cl2 C13 


^c 23 e«, 

C13 


fer 


S12 c 23 

Cl2 c 13 


^23 e« 

Cl3 






C? 3 S 2 3C23 


c 13 c 23 ~ 


f r 2 / 3 (si2Si3C23 e~ i<5 + Ci2S2 3 ) 2 e~ i7 



7"2/3(Sl2.Sl3C23 e ^ + C 12-§23) (^12^13^23 g ^ - Ci 2 C 23 )e * 7 

C13C23 + ^2/3(si 2 si 3 C23 e- i5 + ci 2 s 23 ) 2 e- i7 

i<5 



(14) 



C 13 S 23 + r 2/3( s 12 s 13 s 23 e 



Ci2C2 3 ) 2 e- J7 



C13C23 + ^2/3(Sl2Sl 3 C 2 3 e 16 + c 12 s 2 3) 2 e 17 



with r 2 /3 = m 2 /m 3 , while for the inverted one 

feu — c 23 e* 5 ) 

Sl3 

Jer — — S 23 e , 

S13 

- _ ^2/i(-Si2gi3C 23 e~ ,g + Ci 2 g 23 )(gi 2 gi 3 g23 e~ tg - ci 2 c 23 )e~ t7 
' tT r 2/ i(si 2 si 3 c 23 e" i5 + ci 2 s 23 ) 2 e"^ + (ci 2 si 3 c 23 e" i5 - si 2 s 23 ) 2 
(ci 2 Si 3 c 23 e l& — Si 2 s 23 )(c 12 s 13 s 23 e lS + s 12 c 23 ) 
r 2 /i(si 2 si 3 c 23 e"" 5 + ci 2 s 23 ) 2 e"^ + (ci 2 si 3 c 23 e" i<5 - si 2 s 23 ) 2 ' 
^ _ r 2/ i(gi 2 5i 3 g 23 e~ l6 - ci 2 c 23 ) 2 e~' 7 + (ci 2 gi 3 5 23 e~ l6 + si 2 c 23 ) 2 
^2/i(si2Si3C 2 3 e" i<5 + c 12 s 23 ) 2 e-^ + (ci 2 si 3 c 23 e~ iS - Si 2 s 23 ) 2 ' 

with r 2 /i = m 2 /mi. From the first two equations in Eq. (fl4|) . one can see that f efJj will 
be close to f eT in the limit of #i 3 — > and # 23 — >• 7r/4. This fact turns out to be one of 
the origins of the fi-r symmetry as shown in the next section. On the other hands, f e/x 
and f eT in Eq. ( fl5l) always have an opposite sign. This indicates that the inverted case 
cannot be consistent with the \i-t symmetry 2 . 



III. h-t SYMMETRIC LIMIT AND DEVIATION 



In this section, we investigate the Zee-Babu model by considering the /i-r symmetric 
type of the matrix in Eq. (TTTj) as follows 

f A —B B ^ 
-B C D 
BBC 



\ 



(16) 



/ 



which can be diagonalized by the PMNS matrix in Eq. ( [131) with # 23 = ir/A and #i 3 = 0, 
where A, B, C and D are complex values in general. Note that for the matrix in Eq. (fTTT) 
there are only two possible /i-r symmetric limits: (i) Co 



ijj T 



1 (u 



fe 



(CO, 



co' TT and / e „ = f eT ). However, the former condition 



and (ii) uo TT = 1 and f er - Jeu - ^ TT ^ Jeil 
results in m\ = m 3 = or m 2 = m 3 = 0, which must be largely broken in order to fit the 
experimental data. Thus, we will focus on only the latter one. 



2 There is actually a special case in which f^ T = 0. However, once we force the neutrino mass matrix 
to be /i-r symmetric, the theory suffers from the dangerous lepton flavor violating processes, such as 
r — > /ry. 
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A. Normal mass hierarchy 



In the /i-r symmetric limit, the matrix M Q & in Eq. ffTTT) becomes 



M, 



ah 



\ 



* 
* 



-1 



(17) 



and three mixing angles are given by 



23 = -, 0i3 = O, tan 20i 2 = 



t:/f 



The three eigenvalues are found to be 



Ai 

A 3 



12C12 + S? 2 )(^T + 1] 

(2/*s? 2 + 2v / 2/ eAt s 12 ci 2 + c? 2 )0V + r 



(19) 



where either Ai or A 2 always vanishes. Hence, this limit is only consistent with the normal 
mass hierarchy case. For example, the exact TB mixing is obtained from f efl = 1/2 3 , while 
the central value of the mass ratio, which is 7712/^3 — 0.176, corresponds to Q 



11- 



1.27. 



Moreover, in order to fit all central values in Eq. (JTJ), we need to deviate from the /i-r 
symmetric limit and it can be realized with the following data set: 



/ e/i ~ 0.47 - 0.07e i(0 " 27r) , f er ~ 0.51 + 0.07e l(0 " 2w) , 



0J T 



0.85 + (0.00 - 0.06)e 



i(0 - 2tt) 



(20) 



cj^t- ~ -0.95 + (0.19 - 0.23)e i(0 " 2n) , 



where we have varied 5 and 7 from to 2tt based on Eq. ( 1141) . Although the /i-r 
conditions are no longer exact, they remain as good approximations, i.e., u TT ~ 1 and 
f er ~ f efJL . Therefore, we conclude that there probably exists the /i-r symmetry behind 
the TB pattern of neutrino mixings in the case of the normal mass hierarchy. 



Although f efi = — 1 also implies the exact TB mixing, it leads to a vanishing A2 at the same time. 
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4 
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(2,0,1) 


(2,0,1) 











-1 



TABLE I: A particle content and charge assignments 
B. Inverted mass hierarchy 

As mentioned in the previous subsection, the fi-r condition has to be largely broken 
in the case of the inverted mass hierarchy. For instance, the central values in Eq. (CQ) can 
be obtained from 

f efl ~ 7.49e l (° " f eT ~ -6.94^° " ^\ 

u) TT ~ 2.00 + (0.00 - 1.64)e <(0 " 27r) , (21) 
Gj^ t ~ 1.52 + (0.00 - 0.78)e* (0 " 2n \ 

where we have varied 5 and 7 from to 2n based on Eq. ffl5l) . However, contrary to the 
normal hierarchy case, it is difficult to find out possible remnants of the \x-t symmetry 
from Eq. ( T2T|) . i.e., f efJi 7^ f eT and u TT 7^ 1. This suggests that, in the inverted hierarchy 
case, the TB mixing may just be an accidental result due to the suitable parameter 
tunings. 



IV. FROGGATT-NIELSEN-LIKE Z 5 MODEL 

In the previous section, we have obtained the conditions: f eT = f efl and u TT = 1 to 
derive the /i-r symmetric matrix. The former condition is easy to achieve by imposing 
a permutation symmetry, whereas the latter one may not be because u TT includes not 
only Yukawa couplings but also the charged lepton masses. For instance, if we ignore 
the electron mass, the condition becomes m 2 T y* TT = m 2 y* , and it requires a hierarchy 
between Yukawa couplings rather than a permutation relation. To naturally realize the 
fi-r conditions for the normal mass hierarchy case, we adopt the scheme of the Froggatt- 



Nielsen mechanism 19J] and show a specific model based on an Z5 symmetry. The charge 
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TABLE II: Lepton flavor violating processes and their experimental bounds used in our calcu- 
lation. 

assignments of the particles under the symmetries are summarized in Table [H In this 
model, we introduce a gauge singlet scalar <fi with the charge —1 under Z5 and consider 
the higher dimensional operators. Because of the symmetry, at the leading order, the 
Yukawa matrices F a t, and Y a b in Eq. ([3]) turn out to be 

( VeeX Ve,X 2 VeA* \ 

Y ab = 



ab 



feA 2 fer^ 2 ^ 







\-fer\ 2 -fprX J 



\ 



* * y T A 2 ) 



(22) 



where A =< > /A is the suppression factor of the higher dimensional operators with the 
typical energy scale of the Z 5 symmetry, A. Note that due to the symmetry, the charged 
lepton mass matrix is diagonal up to the leading order. We remark that to simplify our 
discussion, we have assumed that the terms like LHLH<p^ n are strongly suppressed by 
an extremely-high energy scale. It is easy to see that if we assume 

m, 



A 



0.06, f er = f efl , = y T 



(23) 



and m e = 0, we obtain the /i-r symmetric Majorana neutrino mass matrix, given by 



M 



ab 



8^f 2 T \ 2 u 



V 



-2/ e 2 ,A 2 (l 

* 
* 



U), 



-f efM X(l + U^ T ] 
-1 



f e/Jl X(l + U^r) 

-1 



(24) 



/ 



where = y^/y*^ and ^ . . 

In the rest of this section, we discuss several experimental constraints including some 
non-standard lepton flavor violating processes. (See Table |TTJ ) As discussed in Sec. Ill, 
the realistic neutrino mixings requires 



e/i — J er 



ffJ,T 

2A 



(25) 



In this case, the strongest constraint on the f a b couplings comes from the — y cy process. 
Combined with the neutrino mixing data, we get the lower bounds on and y m : 

f^X > 0.008, (26) 

> 0.13, (27) 

and the allowed range for the singly charged scalar mass: 

10 2 GeV < M h < 10 4 GeV (28) 

as discussed in Ref. [l^]. Moreover, since it may be natural to take all y a t, to be the 
same order, we can estimate the branching ratios of £~ — y £^£~£~ mediated by the doubly 
charged scalar. The stringent constraint comes form either r — y 3/i or \i — y 3e, given by 



WabV 

Ml 



< 1Q- 7 GeV^ 2 . (29) 



It is clear that these processes can be accessible in the near future with the TeV scale 
doubly charged scalar if y ab ~ 0(0.1). 



V. SUMMARY 



We have investigated the Zee-Babu model and tried to find out a possible flavor 
symmetry behind the TB neutrino mixing matrix. We have found that there probably 
exists the /z-r symmetry in the normal neutrino mass hierarchy case, but the TB mixing 
may be accidental in the case of the inverted one. We have also attempted to derive 
the fi-T symmetric neutrino mass matrix with a Froggatt-Nielsen-like Z5 symmetry and 
estimated several constraints coming from lepton flavor violating processes. 
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